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Results are presented for the phase separation process of a binary mixture subject to an uniform 
shear flow quenched from a disordered to a homogeneous ordered phase. The kinetics of the process is 
described in the context of the time-dependent Ginzburg-Landau equation with an external velocity 
term. The large-n approximation is used to study the evolution of the model in the presence of a 
stationary flow and in the case of an oscillating shear. 

For stationary flow we show that the structure factor obeys a generalized dynamical scaling. 
The domains grow with different typical lengthscales and R± respectively in the flow direction 
and perpendicularly to it. In the scaling regime R± ~ t"^ and Rx ~ 7t°'°' (with logarithmic 
corrections), 7 being the shear rate, with = 5/4 and qx = 1/4. The excess viscosity Arj after 
reaching a maximum relaxes to zero as 7~^t~"^''^. A77 and other observables exhibit log-time periodic 
oscillations which can be interpreted as due to a growth mechanism where stretching and break-up 
of domains cyclically occur. 

In the case of an oscillating shear a cross-over phenomenon is observed: Initially the evolution is 
characterized by the same growth exponents as for a stationary flow. For longer times the phase 
separating structure cannot align with the oscillating drift and a different regime is entered with an 
isotropic growth and the same exponents of the case without shear. 
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I. INTRODUCTION 

The kinetics of phase separation of a disordered system quenched into a multiphase coexistence region has been 
extensively studied in the last years 0]. The main features of the process are well understood: After an early stage 
during which ordered domains of the equilibrium phases are formed the segregation proceeds in the late stage by 
coarsening of ordered regions according to the power law growth R(t) for the average domains size. In binary 

liquids, the existence of several regimes characterized by different exponents a, due to the presence of various growth 
mechanisms, is well established In these regimes the pair correlation function C(r, t) verifies a dynamical scaling 
law according to which it can be written as C(r, t) ~ f{r/R), where f{x) is a scaling function 

^From the theoretical point of view the most relevant progresses have been achieved in the framework of the 
continuous approach based on the Cahn-Hilliard equation with a Ginzburg-Landau free energy functional, the time 
dependent Ginzburg-Landau (TDGL) model. Within this approach, which neglects hydrodynamics, the properties of 
the phase-separation kinetics can be efficiently studied by means of numerical simulations or analytically in the context 
of approximate theories, among which the so-called large-n limit (one- loop approximation). For a vectorial system 
with an infinite number of components n, indeed, the TDGL model is exactly soluble. The one-loop approximation is 
known to provide a mean-field picture of the phase-separation process which captures the essence of the phenomenon 
at a semi-quantitative level [||. 

In this paper we study the process of phase separation in a binary mixture subject to an uniform shear flow. 
When shear is applied to the system the time evolution is substantially different from that of ordinary spinodal 
decomposition. We consider both a stationary flow and an oscillating shear. 

A stationary flow induces strong deformations of the domains formed after the quench which become 

anisotropic and stretched along the flow direction. Consequently the growth rate along the flow is larger than in 
the other directions. In some experiments a power law increase of the typical size of the domains is observed and 
a value Aa = ax — ct± in the range 0.8 1 for the difference between the exponents in the flow and in the shear 
directions is reported |^,^. Two dimensional molecular dynamic simulations flnd a slightly smaller value [ic| ]. In 
other experimental realizations, when the shear is strong enough, stringiike domains have been observed to extend 
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macroscopically in the direction of the flow 1 11 1 preventing complete phase separation. In general, the scaling behav- 
ior of sheared systems is not clearly understood and the very existence of a scaling regime in different experimental 
systems is questionable. 

In a previous paper |]l2| we have shown that the numerical solution of the one-loop approximation to the TDGL 
model for phase separation under shear exhibits a generalized scaling symmetry characterized by Aa = 1. In the 
scaling regime the structure factor and other observables exhibit the interesting feature of an oscillatory pattern which 
can be related to a mechanism of storing and dissipation of elastic energy where domains are stretched and broken 
cyclically. This new effect has been shown to persist up to the longest available time of our computation and represents 
the hallmark of a complex dynamical pattern induced by the presence of the shear. In a recent paper Rapapa and 
Bray by solving asymptotically the one- loop equations, confirmed analytically the existence of a (multi)-scaling 
symmetry; in the long time limit, however, they do not recover the cyclical pattern described insofar and they infer 
'that the observed oscillations are slowly-decaying preasymptotic transients'. Since their solution is obtained in the 
infinite time limit, then, a reference theory for the description of this remarkable phenomenon is lacking. 

Given that the one-loop approximation is a mean-field solution in spirit the natural question of its accuracy for 
the description of the original model arises. A numerical analysis of the exact TDGL model has been performed 
recently in [Q where it is shown that the global picture of the one-loop approximation is adequate. In particular the 
oscillatory pattern is recovered. The existence of a scaling symmetry and the determination of the related exponents, 
however, has not been clearly established numerically mainly due to finite size effects limitations. The actual value 
of the growth exponents can be inferred by scaling |13| or renormalization group arguments to be a± = 1/3, as 
in the case without shear (we stress the fact that hydrodynamic effects are neglected in this model), and ax = 4/3. 

The shear also induces a peculiar rheological behavior. The break-up of the stretched domains liberates an energy 
which gives rise to an increase Ary of the viscosity . Experiments and simulations show that the excess viscosity 

Ar] reaches a maximum at t = and then relaxes to smaller values. The maximum of the excess viscosity is expected 
to occur at a fixed 'yt and to scale as Ar]{tjn) ~ 7"*^ l^iD- Simple scaling arguments predict = 2/3 ||], but different 
values have been reported [|| . All these features are adequately described by the TDGL already at the one- loop level. 

In this paper we present a complete scenario of the behavior of the TDGL model for phase separation in a shear 
flow in the framework of the largc-n approximation. The behavior of the system is studied along the whole time 
history, from the instant of the quench onward, both in the presence of a steady flow and in the case of an oscillating 
shear where interesting effects are undercovered. Results are presented for two and three dimensional systems. 

This paper is schematically divided as follows: In Sec. 2 we specify the model and introduce the one-loop approxi- 
mation that will be studied thoroughly in the following Sections. Section 3 is devoted to the analysis of the behavior 
of the model subjected to a steady flow. In Sec. 4 the dynamics in the presence of an oscillatory shear is considered. 
In Sec. 5 we present a discussion of the results, debate some open problems and draw our conclusions. 



II. THE MODEL 



The binary mixture is described at equilibrium by a Ginzburg-Landau free-energy 

= Jd'^r^^'+^-^^ + ^lW^l^} (1) 

where if is the order parameter which represents the concentration difference between the two components. The 
values of b, k are positive for any temperature T of the fluid. The parameter a separates stable states of the blend 
with a > ac(T) {ac{T) < 0), from the thermodynamically unstable states with a < adT) where the system phase 
separates. The time evolution of the order parameter is given by the convection-diffusion equation 

^ + v.(^«) = rv2^ + ^ (2) 

where the gaussian stochastic field 77, with expectations 

{ri(f, t)r]{r',t')) = ^2Tr\7^6{f- f')6{t - t') (3) 

describes thermal fluctuations |^. In Eq. (|^) F is a transport coefficient and the symbol (...) indicates the ensemble 
average. The external velocity field here considered is of the form 

V = "fyex (4) 
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where 7 is the spatiaUy homogeneous shear rate that may however depend on time, and is a unit vector in 
the flow direction. In the following we will consider a quench from an uncorrelated isotropic high temperature initial 
condition at the critical composition, i.e. with ((^(r, 0)) = and < <p(r, 0)</3(f^, 0) >= Ai5(r — r*). The main observable 
for the description of the phase-separation kinetics is the structure factor 

C(fc, t) =< Lp(k, t)ip{-k, t) > (5) 

where ip{k, t) is the Fourier transform of the field (/?(r, t) solution of Eq. (|^). In the high temperature initial state we 
consider one has C(fc,0) — A. 

The cubic term in the derivative 5J- /Stp prevents an exact solution of the Eq. (^, as in the case without shear 
However a soluble model is recovered in the one-loop approximation which amounts to the factorization of the cubic 
term of Eq. (^) as 

It is possible to show that the substitution (H) becomes exact in models with a vectorial order parameter when 
the number n of its components becomes infinite. Since {ip^) = S{t) does not depend on space, due to translational 
invariance, the substitution (|^) formally linearizes the theory. The large-n limit is a well developed approximation 
scheme in statistical mechanics which have been applied to different contexts [Q: Its validity and limitations are 
nowadays rather well understood |p^ . 

In the large-n approximation the dynamical equation for C(fc, t) is; 

^^m^ - ^^^^^^ = -^'[fc' + Sit) i]c{lt) + er (7) 

where the function S{t) is self-consistently given by 

S{t)-i J^Mk,t) (8) 



|fc|<9 



and 5 is a high momentum phenomenological cut-off. Notice that in Eq. (Q) the parameters of the free-energy (|l|) and 
the mobility T have been eliminated by a redefinition of the time, space and field scales. The rheological properties 
of the mixture are described in terms of the shear stress 

r dk 

J\k\<q [^■^r 

and of the first and second normal stress differences 

^Ni^ I -—^[kl-kl]C{k,t). (10) 

^|fc|<9 (271-)" 

and 

r dk 

^^-= i T^Akl-kl]C{Kt). (11) 

J\k\<q (^^) 

For vectorial systems with 71 > rf (d is the spatial dimensionality) topological defects are not stable For large n, 
therefore, domains of the equilibrium phases are, strictly speaking, absent. Nevertheless, since from the solution of 
the one-loop equations presented below it is possible to identify characteristic growing lengths Rx (t) and R± (t) in the 
fiow and in the other directions it is natural to interpret these quantities as the trace of the domains size after the 
one-loop approximation procedure has been performed. In the following we will always use the word domains in this 
broad acception. 

III. STEADY SHEAR 

In this section we consider the case of a constant shear rate 7. Eq. (|^) can be formally integrated, yielding 



3 



Jo 

where 

IC{u) ~ k + jkxUCy (13) 

and Cy is the unit vector in the shear direction normal to the flow. For steady flow it is usual to define the excess 
viscosity as 

Arm = (14) 



A. Analytic solution in the short and long time limit 

The consistency condition (^ cannot be worked out along the whole time history of the system. For this reason in 
the following Sections the model Equations will be solved numerically both in d ^ 2 and d — 3. However the model 
can be solved in the short and long time limit. 



Short times 



For short times the linearized theory developed originally by Cahn and Hilliard for the situation with 7 = can 
be extended to the present case. This amounts to neglecting the quartic term in the local part of the free energy (|^) 
since in the initial high temperature state the order parameter is small. With this approximation the solution of 
Eq. (^ reads 







This approach applies to the original model and to the large-n approximation as well because non-linear terms are 
neglected. It is well known that the linear theory describes the very initial transient of the phase-separation process, 
when domains are still forming. In this time domain the behaviour of the system in the presence of the flow is more 
interesting than in the simple case of an immobile fluid. A plot of the structure factor ( p^ is presented for a two- 
dimensional system in Fig. (|l|) for the case 7 = 1 and T — 0. Initially, when domains are forming but the shear flow 
has not yet produced sensible effects, the structure factor evolves assuming the typical structure of a circular volcano, 
similarly to what happens in the case without shear. At "ft c^i 0.5 the anisotropy induced by the shear produces 
a deformation in the proflle of the edge of the volcano from a ring-like geometry into an ellipse, whose major axis 
forms with the positive direction of the ky axis an angle of approximatively 45° (see Fig. (|l)). At the same time 
small dips start to develop in the edge at the ends of the axes of the ellipse and four peaks can be clearly observed 
at 7i ~ 2. As time goes by, the angle formed by the major axis of the ellipse with the fcy-direction decreases and 
the dips in the proflle of C{k,t) along the major axis develop until C{k,t) almost consists of two separated foils, at 
jt ~ 4, when a couple of peaks prevails. The same initial pattern is also observed [ p^ by numerically solving the 
full model equation (^. At later times, however, the presence of the non-linear terms becomes fundamental and the 
linear theory breaks down, as in the case without shear. It is important to stress the fact that the presence of four 
peaks in the structure factor is exhibited already at the linear theory level of approximation. We will see in the 
following sections that the very existence of a multiply peaked C(fc, t) produces a rich dynamical pattern originating 
an oscillatory phenomenon. 



Long times 



The self-consistency condition (^) has been worked out explicitly in the long-time domain in [Q. It is found that 
the model is interested by a multiscaling symmetry, as in the case without shear |2l| | , characterized by the growth of 
the characteristic lengthscales as 
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i^.-7(^)^ (16) 

and 

R,^.(±^)i (17) 

in the directions of the flow and perpendicular to it respectively. The excess viscosity and the normal stress differences 
behave as 

Arjit)r^j'\^ff^ (18) 



ATVi - AA^2 (^)^ (19) 
The same behaviors (apart from logarithmic corrections) is obtained in p2| by means of a scaling ansatz. 



B. Numerical solution 

We present in this Section the results of the numerical integration of the large-n equation (^ which allows to follow 
the whole time history of the phase-separation process. We restrict ourselves to the case with T — 0. An Euler 
first order discretization scheme has been implemented in d = 2 and d = 3 on d-dimensional lattices with 201 mesh 
points per each direction. For long times the structure factor is strongly peaked around typical wavevectors which 
move towards zero as time goes on (see Fig Given that the support of C(fc, t) also shrinks to zero it is possible 
to greatly improve the quality of the numerical computation by using a self-adaptive mesh algorithm that follows the 
evolution of the support of the structure factor. We have solved Eq. (0) for various values of the shear rate 7 in the 
range [10~^, 10~^]. We found that the qualitative behavior is the same for all the values of 7 considered. /^From the 
knowledge of the structure factor we compute the characteristic lengths R{t) as 




where 

J dkC{k, t) 

and the same for the other directions. 



d^2 

The behavior of C(fc, t) is shown in Fig. (H) for 7 — 0.001. Initially the evolution of the structure factor is resemblant 
to the one observed in Fig. (^ where the linear theory for C(fc, t) was plotted. Later on, however, the linear theory 
fails because the non-linearities become relevant, and the long-time regime is entered. This is characterized by the 
shrinking of the support of C(fc, t) towards the origin with different rates for the shear and the flow directions so that 
the tilt angle, namely the direction along which C is aligned, decreases in time. The structure factor is divided into 
two separated foils which are symmetric due to the property C(fc, t) = C{~k, t). In each foil two distinct peaks can be 
observed located at {kxi,ky^) and (^2:2 1^1/2) with \kx^ \ ^2\kx2 \ and |fcj,J ~2\kyr^\. Their heights change in time. The 
first peak to prevail is that located at (kx^ ky-^), while the other peak dominates later. As time elapses the two peaks 
are observed to prevail alternatively. This oscillatory behavior continues up to the longest times of our computations. 

In Fig. (^ the quantities In t)^/'^Rx{t) and {'^hit)^/*Ry{t) are plotted against the strain jt. According to 
Eqs.( ^PT| ) for long times these quantities should collapse, for different values of the shear, on two power- law master- 
curves with exponents 5/4 and 1/4, respectively. Here we observe that the collapse is indeed good, but the predicted 
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power-law behavior is modulated by an oscillatory pattern. These oscillations are observed to be periodic on a log-time 
axis and persist up to the limit of the computational time. 

We now consider the rheological behavior of the mixture by plotting in Fig. the quantity (7/ Int)^/^ A77(t) against 
the strain. This quantity reaches a maximum at 7t ~ 3.5 and then decrease as also found in experiments For long 
times Eq. ( p^ would predict a data collapse for different 7 on a single power-law master-curve with exponent —3/2. 
Here the situation is similar to the previous figure, in that the predicted behavior is modulated by log-time periodic 
oscillations. On the bases of simple scaling arguments the maximum of the excess viscosity Ar](tm) is expected to 
occur at a fixed -ft and to scale as Ari{tm) ~ 7"", with v = 2/3 These arguments do not directly apply to the 

one-loop approximation since, due to the mean field nature, the exponents are different. The asymptotic solution ( p^ ) 
is not adequate to this early stage effect. The 7 dependence of Ari{tm) is plotted in the inset of Fig. ^ showing that 
a power-law behaviour with i' ~ 0.6 is obeyed, in partial agreement with the aforementioned scaling arguments. 

In Fig. (|^) we report the numerical results for the first normal stress by plotting (7lni)~^/^A7Vi against -ft with 
7 — 0.01. We find that AiVi(i) scales asymptotically as predicted by Eq. ( pj| ) again modulated by an oscillatory 
pattern. 

The periodic oscillations observed in all the physical observables are due to the competition between the different 
peaks of C{k,t). Let's refer to the behavior of the excess viscosity to understand how this competition affects the 
rheological quantities, using the features of the structure factor to obtain information about the domains evolution 
under the action of shear. Arj reaches its first maximum when the shape of C{k,t) is such that the peak located 
at (kj.-^, ky-^) prevails and the difference between the height of the two peaks is maximal. At this time the domains 
are elongated by the flow and there is a prevalence of thin domains in the system. As these string-like domains are 
stretched further, they eventually break up into two or more domains, dissipating the stored energy. This has two 
effects: the excess viscosity decreases and, on the other hand, the thick domains, which have not yet been broken, 
prevail. In this situation the other peak of C'{k, t) (which is located at (fca;^, ky^) and represents the smaller features) 
grows faster until it prevails and A77 reaches a minimum. This behavior is reproduced with a characteristic frequency 
in log-time. Recently, a similar behavior has been observed in the numerical simulation of the full model Equation (0) 




In this section we report the result of the numerical solution of Eq. in d = 3. In Fig. (|^) the time evolution of 
the structure factor in the special planes fc^, = 0, /cj, = and fc^ = is shown for 7 = 0.001. In the plane k^ — C{k, t) 
behaves analogously to the previously discussed two dimensional case. The structure factor on the plane ky = gives 
information relative to the observation of the system along the shear direction: No velocity gradient is present in 
the plane perpendicular to this orientation, but there are different velocities in the x and z directions. This allows 
to explain the observed behavior which is rather different from the one observed at k^ = 0. The structure factor 
develops initially a circular volcano, as without shear. The edge of the volcano is progressively deformed by the shear 
into an ellipse with the major axis along the kz direction. The dips in the edge of the volcano at values of fc^, ~ 
develop with time so that at -ft ~ 1, C{k, t) is made of two foils but these are not completely separated. During the 
time evolution the axes of the ellipse shrink; the decrease is faster along the k^ direction. The two foils are never 
completely separated and the angle formed with the k^ direction is zero, as observed in experiments At 7t ~ 5 
two well-formed peaks start to develop and grow on each foil of C{k,t). These four peaks have the same height and 
their relative heights do not change in time, as it can be seen at -ft = 20 in the picture, differently from the situation 
on the kz = plane. 

In the kx = plane the shear has no effect at all and the structure factor remains circular during its evolution. 
The computed behavior of Rx{t) and Ry{t) is similar to that of the two-dimensional case. We also find Rz{t) ^ Ry{t), 
as expected. 

We report in Fig. ^ the plots of {-f/lnty/'^Ar]{t), {-flnt)-^^^ANi{t) and -{-flnt)-^/'^ AN2{t) as functions of 7i. 
It appears that the rheological quantities still have amplitudes which are modulated by log-time oscillations which 
are in phase among them. The origin of such oscillations has to be found again in the oscillations of the peaks of the 
structure factor in the plane kz = 0- Since the support of C(A:, t) shrinks towards the origin faster in the kz than in the 
ky direction the second normal stress difference AN2{t) is negative. This is in accordance with general experimental 
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IV. OSCILLATING SHEAR 



In this Section we consider the case of a time-dependent shear rate with 



j{t) = 7o cos Lot 



(22) 



This situation is of great experimental relevance expecially for probing the viscoelastic properties of the phase sepa- 
rating binary mixture. 

We solved the Eq. (0) numerically in d = 2 using the same numerical scheme as in the case of steady shear, for 
different values of 70 and w. We will describe below the case 70 = 10""^, r = 27r/w = 6 x 10^. The time evolution 
of the structure factor in the first cycle of 7(t) is shown in Fig. (^). The dynamical pattern is analogous to the one 
with 7 = const, for times t < t/4 as it can be seen at "fot = 1.5. Then the time dependent velocity field modifies 
the behaviour of the blend with respect to the case of a steady flow. In particular, at the end of the first oscillation, 
the four peaks of C{k,t) are located at comparable distances from the origin of the fc-space differently from what 
observed in Fig. (||) at jt — 6. The two highest maxima at 70^ = 6 in Fig. (H) are characterized by \ky\ ^ \kx\- 
During the later time evolution these peaks grow and move towards the origin. The position in the fc-plane of the 
other peaks rotates back and forth cyclically along an approximatively circular path. The radius of this trajectory 
shrinks towards the origin at a rate comparable with that of the position of the other peaks. In the asymptotic regime 
the four peaks have approximatively the same height and the cyclical rotation of the peak position persists. This is 
shown in Fig. (^) where the configurations of the structure factor are shown at each quarter of oscillation of the shear 
rate in the asymptotic stage. 

In Fig. ( p^ the evolution of the characteristic lengths Rx{t), Ry{t) is plotted against 70^. We also plot, in the inset, 
the time average of these quantities over a period r, in order to smoothen out the superimposed oscillations. Here we 
observe, for times t < t, growth laws analogous to the steady shear case, namely Rx{t) ~ t^^^ and Ry{t) ~ t^/"*. The 
growth exponent of Rx changes smoothly, from i ~ t onward, from 5/4 to the asymptotic value 1/4 which is reached 
at 7oi ~ 80 when all the four peaks of the structure factor have the same height. The gradual crossover of ax from 
5/4 to 1/4 can be better observed for larger values of r, since the regime with ax — 5/4 persists for a longer time. 
This is shown in Fig. (|ll|), where the evolution of Rx{t), Ry{t) is plotted against 7oi for r = 5 x 10'"'. For small r, 
instead, Rx and Ry grow with the same exponent 1/4 from the beginning. 

These observations suggest the following physical interpretation: for t < r/2, since 7 does not change sign, the 
evolution of the blend is comparable to the case with a constant shear rate. In particular, if t is sufficiently large to 
exceed the initial stage when domains are forming, the power growth laws described in Sec. 3 for Rx, R± are observed 
with ax = 5/4 and a± = 1/4. On timescales much longer then r, however, the network of the larger domains cannot 
be efhciently tilted along the flow orientation which changes periodically its sign. This is confirmed by the behaviour 
of the two peaks with \ky \ ^ \kx\ whose position in the fc-plane moves toward the origin but does not cross the fc^, = 
plane, as it would be the case if the orientation of the domains corresponding to these peaks were reversed. In this 
situation the difference Aa = 1 between the exponents in the flow and shear direction cannot be sustained, because 
the larger domains are not directed along the flow orientation at all times, and a growth law with the same exponent 
1/4 in all the directions is obeyed. It is interesting to notice that the other peaks, which represent smaller domains 
formed by the break-up of the larger ones, crosses the ky — plane during their rotation every half period of 7. This 
suggests that these features are tilted by the oscillating shear and follow the flow orientation. Then we expect to 
observe in a real blend two type of domains which respond differently to the oscillations of the flow: a network of large 
and elongated structures which maintains the orientation imposed during the first half period of 7 and a multitude 
of more isotropic features, generated by the break-up of strained regions, which oscillate following the flow. 

For studying rheological properties it is customary to introduce a complex viscosity rj* = rj' ~ irf which is 
related to the shear stress by 



<^xy(t) = Joiv' cos Lut + rj" siriLut) 



(23) 



when Eq. (22) holds. It is also useful to consider [g4| another representation of the shear stress given by 



'^xyit) — C sin(a;i -I- (f)) 



(24) 



The connection between Eqs. (|23| ) and (|4|) is given by 



C = 7oV'?'^ 



(25) 



and 
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tan0=^. (26) 
rj" 

By defining 7*(t) = 7oe*"*, we can write Eq. ( p3| ) as 

a,,(<)=i?e[77*7*(t)] (27) 
In order to relate the real and imaginary parts of the viscosity to physical quantities Eq. ( p3| ) can be casted as 

<rxy{t) = -niit) + G I i{t')dt' (28) 

where r] — rj' , G — uii]" and the identity sina;t = lu cos tot' dt' has been used. 

The coefiicient rj in the r.h.s. of the Eq. ( p8| ) multiplies the portion of the shear stress in phase with the shear rate 
and represents the viscosity of a viscoelastic fluid. The integral of the second term of the r.h.s. of the Eq. ( p8| ) can 
be identified with the shear strain present in the mixture at time t. The coefficient G is therefore the effective elastic 
shear modulus of the fluid. Pure viscous behavior corresponds to G = ((/) = 7r/2), pure elastic behavior to 77 = 
(0 = 0) H. 

In order to compute rj and G during the phase separation we calculated by numerical integration the shear stress 
using its general definition (j^). By writing 

Cxyit) — A cos ujt + B smujt (29) 

it follows that -q — A{t)/^o and G ~ ujB{t)/jo- In general, A and B depend on time. During a single shear oscillation, 
however, we expect that the Eq. ( |29| ) holds as a good approximation with constant values for A and B. In this way 
<^xy{t) is expressed in terms of the first two coefficients in a Fourier series expansion over the interval of scaled time of 
duration 27r. The values we obtain may be referred to the time that locates the middle of the interval. Thus we get 



1 \ \ 1 



2 



in-2-K 



?y m-- r = / axy{t/uj)cost dt (30) 



Tol" J(m-l)-27r 



1 \ \ />m-27r 
1 \ \ LO 



G ( ( m — - I T j = j axy{t/uj) suit dt (31) 



To"" J(m~l)-2-K 



where m = 1, 2, ... 

In Fig. ( [1^ ) we report the plots of r/ and G against 7oi. The viscosity shows a crossover between a power law decay 
with exponent —3/2 at short times and an asymptotic behavior whose exponent is —1/2. This can be explained 
observing that, for the steady shear case, the dynamic viscosity r/ coincides with the excess viscosity which scales with 
the inverse of the domains volume V . When an oscillatory shear is applied V crosses over from an initial power law 
increase V ~ i"^/^, similar to the one for the case with steady shear (see Fig. (||)), to a slower growth V ^ i^/^, as 
already discussed above for Rx{t), producing a corresponding crossover in 77. 

^From the computed values of r/ and G we estimated the phase angle (f>, which, according to Eqs. (|6|) and (|2^), 

is given hy (j) — arctan (^'q ^^8- ( P^ ) we report the time evolution of (/i as a function of jot- It can be seen 

that (p decreases with time to reach an asymptotic value which is approximately 0.016. Accordingly, the system we 
are investigating shows in the asymptotic stage a behavior which is essentially elastic. The experimental data of ]2^ ] 
confirm this behavior. 



V. SUMMARY AND DISCUSSION 



In this paper we have studied the kinetics of a phase-separating binary fluid, in the presence of a shear flow, by 
means of the TDGL model. It is nowadays well established that the corresponding model with 7 = accurately 
describes the main features of the segregation process in binary alloys, where hydrodynamics can be neglected. In 
viscous fluids, such as polymeric blends, the validity of the present approach is limited to the early stage of spinodal 
decomposition; for longer times one should consider the full hydrodynamic description 

When the shear is applied to a fluid the behavior of the system is profoundly changed under many respects and a 
general agreement on the predictivness of the proposed models is matter of general debate. A discussion on possible 
effects of hydrodynamics is presented in Moreover, the numerical solution of the TDGL model with shear poses 
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serious problems due to discretization limitations and finite size effects and, although some progresses have been 
recently achieved [jl^ , a satisfactory description is nowadays not available. In this scenario it is important to devise 
a simple analytical scheme providing the fundamental tools for the comprehension of the fluid dynamics. A natural 
choice in the field of growth kinetics is the large-n approximation, that has been thoroughly studied in the case 
without shear, where it has proven to give a reliable description of the segregation process, although at a mean field 
level. 

In this paper the behavior of the TDGL model in the one-loop approximation is studied in detail, and the whole 
time evolution of the blend is considered, from the quenching instant onward; the cases of a stationary flow and of 
an oscillating shear have been examined. In doing so we undercover a very rich dynamical pattern, where not only 
some experimental findings are reproduced, but new predictions are allowed. After an early stage, which is accurately 
described by the linear theory d la Cahn-Hilliard, the presence of the velocity field produces an anisotropic power-law 
growth of the characteristic lengths Rx, R± respectively in the flow direction and perpendicularly to it. The value of 
the exponent a± = 1/4 in the directions perpendicular to the flow is the same as in models with vectorial conserved 
order parameter without shear; although the actual value of this exponent is not expected to be accurate for real fluids 
(since, even without shear, the exponent obtained at the same level of approximation is known to corresponds to the 
Lifshitz-Slyozov exponent a = 1/3 for scalar fields) a growth exponent a± unaffected by the presence of shear has 
been obtained also by scaling |]l3| and renormalization group arguments applied to the full model equations and is 
also measured in experiments pT Moreover a difference Aa = 1 between the flow and shear exponents is also expected 
to be obtained by releasing the present approximation |]l^,Q and is observed in some experiments. In the case of a 
stationary flow the anisotropic growth governed by these exponents is observed from the onset of the scaling regime 
onwards. The power law behavior of any observable is decorated by log-time periodic oscillations. These oscillations 
characterize the scaling regime up to the longest simulated time but they are not observed in the asymptotic solution 
presented in [ p^ . Given that log-time periodicity appears to be a rather common feature being observed, besides 
segregating fluids, during fracturing of heterogeneous solids [ p6| , p7[ and in stock market indices |^ for instance, it 
would then be interesting to devise an analytical approach to enlighten the origin of this new phenomenon at least in 
the present model. 

In experiments with real fluid systems carried out by Laufer et al. |^ and, successively, by Mani et al. and 
Migler et al. ||3ll] a double overshoot in the time behavior of the viscosity and of the normal stress is observed and an 
interpretation in terms of break-up and recombination of the domains network is proposed. On the bases of our results 
it is plausible that this double overshoot represents the first part of a log-time periodic phenomenon which could be 
hopefully detected with a suitable experimental setup. In the model we have studied the oscillatory behaviour is due 
to the competition between the different maxima of a four-fold peaked structure factor. The presence of these maxima 
is interpreted is Sec. Ill B as due to the existence of different types of domains and the recurrent prevalence of each 
peak is suggested to be caused by the interplay between these kind of regions. A structure factor with four maxima 
has also been observed in polymer mixtures pl| ; however to our knowledge the connection between the alternative 
dominance of the peaks of C(fc, t) and the overshoots observed in the viscosity and in the stresses has never been 
discussed insofar, perhaps due to insufficient resolution, although an experimental confirmation of this hypotheses 
would be desirable. 

When an oscillating shear is present, the anisotropic regime discussed so far for the steady shear case crosses over 
to an isotropic growth when domains are fully developed. In this late stage, from the analysis of the behavior of 
the structure factor, we conjecture again the existence of two types of domains responding differently to the fiow 
oscillations: the network of elongated structures keeps the orientation assumed during its formation in the early stage 
while small features generated by scission of strained parts oscillate in phase with the flow. In this late stage the 
growth kinetics is regulated by the same exponents as without flow. We are not aware of experiments reporting these 
features: It would be interesting to devise an experimental set-up for testing this prediction. 
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CAPTIONS 



FIG. 1. The evolution of the structure factor in the hnear approximation for 7 = 1. The range of kx and ky varies from -2 
to 2. At 7t = 4 the highest peaks are located at {kx, ky) ~ (0.0, ±0.75), the other two at {kx, ky) ~ (=p0.38, ±0.80). 

FIG. 2. The evolution of the structure factor from the numerical solution of Eq. in d = 2 for 7 = 0.001 and T — 0. The 
range of kx and ky varies, at increasing times, as —0.6 < kx,ky < 0.6 at — 0.05, 1; —0.15 < kx < 0.15 and —0.6 < ky < 0.6 
at 7t = 6, 8; -0.075 < kx < 0.075 and -0.6 < ky < 0.6 at yt = 10; -0.01875 < kx < 0.01875 and -0.3 < ky < 0.3 at -yt = 45. 



FIG. 3. Data collapse (scaling plot) for the domains radii in d = 2. The quantities {-y In t)^^^Rx(t) and ('j In t)^^^Ry{t) are 
plotted against the strain for different values of the shear rate: (A) 7 = 0.0001, (o) 7 — 0.001, (•) 7 — 0.01. The two straight 
hnes have slope 5/4 and 1/4. 



FIG. 4. Data collapse (scaling plot) for the excess viscosity in d = 2. The quantity {'y / Int)^^^ Ari(t) is plotted against the 
strain -yt for different values of the shear rate: (A) 7 = 0.0001, (o) 7 = 0.001, (★) 7 = 0.005, (•) 7 = 0.01. The straight line 
has slope -3/2. The inset shows the maxima of Ari as a function of 7. The slope of the straight line is 0.6. 



FIG. 5. The first normal stress AA'^i multiplied by {'ylnt) ^''^ as a function of -yt. The value of 7 is 0.01. The slope of the 
straight line is -1/2. 

FIG. 6. The structure factor from the numerical solution of Eq. (^ in d = 3 at consecutive times for 7 — 0.001 and T = 0. 
From up to bottom the sections kx — 0, ky — and fez = of the structure factor are shown. The range of kx, ky and kz varies 
as -0.6 < kx, ky, k^ < 0.6 at 7* = 1; -0.0375 < kx < 0.0375, -0.3 < ky < 0.3 and -0.3 < k^ < 0.3 at yt = 20. 

FIG. 7. (7/lnf)^/^Ar7(t) (*), (7 In t)"^/^ AiVi (t) (o) and -(-ylnt)'^^'^ AN2{t) (•) as functions of the strain -yt. The straight 
lines have slope —3/2 and —1/2. 



FIG. 8. The structure factor at each quarter of the first 7(t) oscillation. The range of kx and ky varies, at increasing times, 
as -0.6 < kx, ky < 0.6 at yt = 1.5; -0.3 < kx < 0.3 and -0.6 < fc,; < 0.6 at jt = 3, 4.5; -0.3 <kx,ky < 0.3 at 'yt = 6. 

FIG. 9. The structure factor at each quarter of a single oscillation in the asymptotic stage. The range of kx and ky is 
—0.15 < kx,ky < 0.15 at each time. 

FIG. 10. The size of the domains in the x (upper curve) and y (lower curve) directions are plotted against 7ot for r = 6 x 10'^. 
In the inset the same quantities (* and • for the x and y direction, respectively) averaged over every single period of oscillation 
are plotted at times (m — l/2)7or, with m = 1, 2, .., against 'yot. The two straight lines have slope 1/4. 

FIG. 11. The size of the domains in the x (upper curve) and y (lower curve) directions in the case r = 5 x 10^. The two 
straight lines have slope 5/4 and 1/4. 

FIG. 12. The viscosity 77 (*) and the elastic shear modulus G (•) for r = 6 x 10^ at times (m — l/2)7or with m = 1, 2, ... 
The straight lines have slope —3/2 and —1/2. 



FIG. 13. The time evolution of the phase angle 4> for r = 6 x 10^ at times (m — l/2)7or with m = 1, 2, ... 
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